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Abstract. This paper is a further step in the investigation of the morphology of the color-magnitude diagram of Galactic glob- 
ular clusters, and the fine-tuning of theoretical models, made possible by the recent observational efforts to build homogeneous 
photometric databases. In particular, we ex amine here the calibration of the morphological parameter Whb vs. metallicity, 



originally proposed by Brocato et al. (1998; B98), which essentially measures the color position of the red-giant branch. We 
show that the parameter can be used to have a first-order estimate of the cluster metallicity, since the dispersion around the mean 
trend with [Fe/H] is compatible with the measurement errors. The tight WiiB-CFe/H] relation is then used to show that variations 
in helium content or age do not affect the parameter, whereas it is strongly influenced by the mixing-length parameter a (as 
expected). This fact allows us, for the first time, to state that there is no trend of a with the metal content of a cluster. A thor- 
ough examination of the interrelated questions of the a-elements enhancement and the color- T c s transformations, highlights 
that there is an urgent need for an independent assessment of which of the two presently accepted metallicity scales is the true 
indicator of a cluster's iron content. Whatever scenario is adopted, it also appears that a deep revision of the V — /-temperature 
relations is needed. 

Key words. Stars: evolution - (Stars:) Hertzsprung-Russel (HR) and C-M diagrams - Stars: horizontal-branch - Stars: 
Population II - (Galaxy:) globular clusters: general - Galaxy: halo 



1. Introduction 

The comparison of the observed color-magnitude diagrams 
(CMD) of a star cluster with the theoretical isochrones is the 
best tool we have to tune some fundamental parameters of the 
stellar evolutionary models. Only when we are sure that the in- 
put parameters and the input physics are correct, we can use the 
model to infer some properties of the stellar population (like 
age, helium content, etc.), not otherwise empirically measur- 
able. Clearly, this is a complex job: on one side we want to use 
the cluster stars to tune the models, and on the other we need to 
use the models to infer properties of the cluster stellar popula- 
tion. Not surprisingly, any stellar models must adopt a number 
of assumptions, often not directly supported by observational 
evidence. One of the most uncertain parameters, strongly af- 
fecting the theoretical location of some branches of the CMD, 
like the red giant branch, is the mixing length parameter a. This 
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parameter, in the framework of the mixing-length theory (MLT; 
Bohm-Vi tense 1958 ), determines the efficiency of energy trans- 
port by convection in the outermost layers of a star. For a given 
stellar luminosity, it also determines the exact radius of the star, 
and hence its effective temperature and colors. 

It is well known that, in order to reproduce the typical tem- 
peratures and colors of red giant stars, a is required to have 
a value between 1.5 and 2.0. Also, a value close to 1.7 is re- 
quired for reproducing the solar radius in non-diffusive solar 
models, whereas about 1 .9 is favored when diffusion is taken 
into account. The fact that a is similar for red giants and the 
Sun (VandenBerg et al. [2000| ; Alonso et al. |2000| ) leads to the 
usual approach of calibrating a by means of a solar model - 
i.e. a model with 1 Mq and solar composition, that is required 
to have the solar luminosity and radius at an age of ~ 4.5 Gyr. 
The same a value is then used to model all stars, including 
red giants. However, there is no theoretical justification that the 
same value of a should apply for any star. 
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Fig. 1. M80 (filled squares) vs. M12 (open squares). An offset 
S(V - I) = -0.03 and SV = +1.50 was applied to the M12 
diagram before overlapping it to that of M80 



In this paper, we have attempted to investigate which is 
the best value of a for low mass (globular cluster) stars and, 
overall, whether there is any dependence of a with the cluster 
metallicity. Our results cannot be considered definitive; how- 
ever this paper shows a possible approach to the problem, and 
enlightens all the uncertainties associated to the calibration 
of the mixing length parameter. Also, this work is limited to 
the context of the MLT, which is still a fundamental approach 
adopted in most computations of stellar models. However, one 
should keep in mind that the MLT is admittedly a very approx- 
imate theory, and that alternative approaches (e.g. Canuto & 
Mazzitelli |1991| ; Spruit |1997| ; Ludwig et al. |1999| ) have been 
suggested. 

This project has been stimulated by an investigation car- 
ried out by some of us a few years ago. In the course of a 
photometric st udy o f the Galactic globular cluster (GGC) M80 
(Brocato et al. 1998 hereafter B98), we compared the morpho- 
logical characteristics of its color-magnitude diagram (CMD) 
with those of an ensemble of other globulars for which a CMD 
was available in the literature. That investigation left an open 
question concerning the relative position of the RGB with re- 
spect to the HB as a function of the metallicity. More specifi- 
cally, when the CMD of M80 was compared with that of GGCs 
with similar metallicity, we found a group of clusters (M3, 
M13, NGC 7492) whose members had CMDs overlapping that 
of M80, and another (M12, NGC 1904, NGC 5897) for which 
significant discrepancies were seen. The discrepancies were in 
the sense that while the HBs could be overlapped in a satisfac- 
tory way, the RGB fiducial lines showed a dispersion in color, 
M12 having the reddest branch (cf Fig. 6 and 7 in B98), as ex- 
pected if M12 would be more metal rich than indicated in the 
literature. We also excluded that the discrepancy could be due 
to an age difference. 



In order to have a more quantitative comparison among a 
larger sample of clusters, with different metal content, we de- 
vised a new morphological parameter with the aim of quanti- 
tatively measuring the distance in color of the RGB with re- 
spect to some fixed point on the HB. We selected as a refer- 
ence point the so called "HB turn-down"(HB-TD), since it has 
been demonstrated that the location of this point in the CMD is 
largely unaffected by the cluster age, metallicity, and primor- 
dial He content (cf. the detailed discussion in B98). We fixed 
the HB-TD at (B - V) = and measured the RGB color 
at 0.5 magnitudes brighter than the HB-TD level. As the posi- 
tion of the HB- TD is fixed, Whb just shows the displacement 
of the RGB as a function of the metal content (cf. Fig. 10 of 
B98). We found that the trend with the metallicity is in the di- 
rection expected from the theoretical models, though with a 
large dispersion, particularly evident at intermediate metallic- 
ities (—1.85 <[Fe/H]< —1.50), larger than expected on the 
basis of the measurement errors. 

We were not able to decide whether such an evidence might 
simply be the consequence of the errors in the photometric cal- 
ibration, or a peculiar distribution in the global metallicity of 
the clusters, and deferred further discussion until direct mea- 
surements of a elements and/or a database of CMDs with ho- 
mogeneous calibration would have been available. 

Such database is now available, thanks to the efforts of 
a number of the original investigators of B98, and we are 
now able to re-examine the question, taking also advantage of 
the new theoretical calculations that we have specifically per- 
formed for this project. 

The paper is organized as follows. The datasets are pre- 
sented in Sect. ||. The measurement procedures and the correc- 
tions for differential reddening are described in Sect. 2.1 and 



respectively. Sect. 2.4 re-examines the original question of 



the trend and dispersion of Whb as a function of metallicity. 

Sect. [f] deals with the comparison of the VandenBerg et 
al. ( |200q ; V00) and the Girardi et al. ( |2000fc GOO) theoretical 
isochrones to the data. Due to the lower degree of sampling 
of the GOO isochrones, a more detailed description of the turn- 



down identification and its measurement is offered in Sect. 3.2 



The keys to the interpretation of the observed vs. computed 
trend of Whb are offered in Sect. In particular, the depen- 
dence on metallicity, mixing length, age, helium content, a- 
enhancement, and the T ff-color transformations are examined 
in Sects. patp^ptpjpl^ndPl 

Whether our data can be used to constrain the mixing length 
parameter is examined in Sect. |5| Furthermore, we suggest in 
Sect. ^ that the parameter can be employed as a metallicity in- 
dicator, with accuracies comparable to other more widely used 
photometric indices. Our summary and conclusions are given 
in Sect. f| 

The nomenclature deserves a final note. We will be using 
several symbols for our parameter throughout this paper. When 
talking about the parameter in a general way, we will use the 
old symbol Whb, but when referring to measurements specifi- 
cally made for the (B — V) or (V — I) colors, we will use the 
symbols W^ v and , respectively. 
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2. The data set 

For this paper we used two new, photometrically homogeneous 
databases that we have recently created: 

- The V, I band "ground-based" dataset collected by 



Rosenberg et al (2000a, 2000b ) used for the GGC relative 



age project by Rosenberg et al. (1999); 



- the "HST snapshot" dataset (Piotto et al. 2002), based on 



F439W (B) and F555W (V) WFPC2 images of the core 
of all the GGCs with (m - M) B < 18.0, which has al- 
ready been used to co nstrain a number of param eters in the 
mode ls (Zoccali et al. |2000j Piotto et al. |2000t Bono et al. 
2001b. 



Both databases are available to the community on 
the Web pages of the Padova Globular Cluster Group at 
http://menhir.pd.astro.it. 

For the present project, we used 26 clusters in the B, V 
bands from the HST snapshot database, spanning a metallic- 
ity interval —2.2 <[Fe/H]< —0.5, and 31 clusters from the 
ground-based dataset, in the V, I bands, and covering a metal- 
licity interval -2.3 <[Fe/H]< -1.1. 

The most important property of these databases for the aim 
of the present project is thei r ph otometr ic internal homo geneit y 
(cf. Rosenberg et al. |2000H and |2000bj and Piotto et al. |2002| ). 
This homogeneity removes the main uncertainty in the com- 
parison done by B98. And, indeed, as exemplified in Fig. |lj the 
two new CMDs of M12 and M80 perfectly overlap, giving the 
first confirmation that the spread noted by B98 was mainly due 
to the photometric in-homogeneity of the data collected from 
the literature. 

2. 1. Measurement of the Wub parameter 

In order to compute the value of Whb, we first measured the 
position of the HB turn-down for each cluster, and then the 
color of the RGB at 0.5 magnitudes above the HB-TD. 

The HB-TD position is, using the definition of B98, the 
place where the true color of the HB is zero. In order to find the 
TD on the observed diagrams, we first defined a fiducial HB, 
that extends from the faint-blue tail to the bright-red clump. For 
the HST dataset, NGC 2808 was used, since it al ready satisfies 
our requirements (see e.g. Fig. 4 of Bedin et al. 2000 ). On the 
other hand, no cluster in the ground-based sample has such an 
extended HB, so we constructed a mean fiducial HB, following 
the procedure adopted in Rosenberg et al. (1999 1. Briefly, we 
started with the HB of NGC 1851, which has a bimodal HB, 
and then extended it to the blue or the red by comparison with, 
respectively, more metal poor and more metal rich clusters. 

The theoretical HBs from VandenBerg et al. (2000; see 
Fig. ||) were then fitted to the fiducial HBs, and the TD-HB 
fixed at (B - V)=0.0, and (V - i>0.0, as in B98. For each 
cluster, we applied a color and a magnitude offset until a sat- 
isfying visual superposition to the reference HB was obtained 
(cf. Figs. ||, Q). Finally the TD color and magnitude for each 
cluster are those of the mean fiducial plus the two offsets. A 
special case is that of NGC 6388 and NGC 6441, whose HBs 
are evidently sloped, such that the blue portion just before the 




Fig. 2. The TD position is shown on the V00 theoretical HBs in 
the case of the (B — V) color (upper panel) and (V — I) color 
(lower panel). Each line corresponds to the theoretical ZAHB 
sequence for a given metallicity. Each one of them has been 
vertically shifted in this plot, so that they all coincide at the TD 
point. In the two panels, the empirical fiducial HBs have been 
represented as well. They can be recognized as the solid curves 
ending at colors w 0.7, and slightly fainter than the theoretical 
loci. 



HB-TD is brighter than the red portion (see Piotto et al. 2002 ). 
There are enough stars on the blue side, around the HB-TD, to 
allow a determination of the HB-TD position, without relying 
on the red part of the HB. However, if we tried to reach an over- 
all agreement from blue to red, then a slightly fainter HB-TD 
would have been found. In turn, this implies a smaller Whb> 
and one can see in Fig. [7] and the following, that the represen- 
tative points of the two cluster would move toward the general 
trend defined by the isochrones. 

In order to measure the color of the RGB, we first obtained 
by hand an approximate color of the RGB region 0.5 magnitude 
brighter than the TD level. The final color was computed as the 
median color of all the stars comprised in a rectangular region 
around this first RGB position estimate. The box dimensions 
were slightly varied from cluster to cluster, so that the whole 
RGB color extension was comprised. The "vertical" size was 
chosen as to ensure that a statistically significant number of 
stars entered the computation. Typical values are 0.06 magni- 
tudes in color and 0.12 magnitudes in V. The Figs. |] and Q 
show the final step of the procedure for NGC 5694 (HST data) 
and NGC 6656 (groundbased data). 

2.2. Correction for differential reddening 

Since Whb represents a difference in color, and since the two 
points are about 1 magnitude apart, we cannot neglect the ef- 
fect of the wavelength dependence of the reddening. It is known 
that, given a mean reddening Eb-v, the value of the absorption 
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Fig. 3. An example of measurement of W^ v , for NGC 5694. 
This figure shows the mean fiducial HB in (B — V), V (large 
open circles), the position of the TD, and the corresponding 
color of the point on the RGB which is 0.5 magnitudes brighter 
than the TD. The inset shows the histogram of the color distri- 
bution on the RGB, and the solid curve is a Gaussian computed 
assuming the calculated mean and dispersion of the data. 
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Fig. 5. All the data for W^ v , as a function of metallicity, 
adopting the CG97 (top panel) or ZW84 (bottom panel) metal- 
licity scales. Filled triangles identify clusters w hose m etallicity 
on the ZW84 scale has been taken from Harris (1996). 




Fig. 4. Same as Fig. ||, showing an example of measurement of 



Whb 7 for NGC 6656 - 



only as a function of the temperature, but also on star's gravity 
and metallicity. The results of GR95 are consistent with Olson 
dl975|). 



We used the results from both studies to correct the Whb 
parameter measured in the B, V bands. In order to correct the 
Whs parameter obtained in the V vs. (V — I) CMD, we used 
the GR95 tables, thus ensuring the homogeneity of the results. 
More details on our calculations are given in Appendix ^| 

2.3. Metallicity scales 

Most values of the iron abundance for the single clusters have 
been taken from Rutledge et al. ( [1997} RHS97). More explic- 
itly, we have taken the [Fe/H] values that RHS97 obtained by 
calibrating their Call triplet equivalent widths onto either the 



Zinn & West (|1984f ZW84) or the Carretta & Gratton ( |1997 
CG97) metallicity scales. For a few clusters, there are no esti- 
mates of the metallicity provided by RHS97, and in such cases 



the [Fe/H] values on the ZW84 scale are from Harris (1996), 
and from Carretta (private communication) for the CG97 scale. 
We have attached a typical error of 0.15 dex to the metallicities 
in this class. 



Ay = Rv x Eb-v varies with the star color (temperature). In 
the case of heavily reddened clusters, we must also take into ac- 
count the dependence of Ry on the absolute reddening Eb-v- 
A first in vestigation on these effects was carried out by 
Olson ( 1975 ), who propos ed a co rrection for the (B — V) color. 
Later, Grebel and Roberts ( |l995j GR95) extended this study by 
determining the trends of R, Ay, colors, and reddenings, not 



2.4. The Web parameter vs. metallicity 

As in B98, Figs. | and |display the W^ B V and 1 param 



'HB a " u yr HB 

eters as a function of [Fe/H], both on the Zi nn an d West ( |1984 , 
lower panel) and the Carretta and Gratton ( 1997 , upper panel) 
metallicity scales. The same data is presented in Table 1. The 
first remarkable result, at variance with B98, is the low disper- 
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Fig. 6. The same as Fig. || but for W^ 1 . 



sion of the data points, fully compatible with the error bars. 
The dispersion at intermediate metallicities noted in B98, is no 
longer present. As originally suspected, the anomalous trend 
in the B98 data was likely due to calibration errors present in 
their CMD database, which was a simple collection of litera- 
ture data. Again, this result shows the importance of using a 
photometric homogeneous database in deriving the properties 
of the stellar population of star clusters. 

The linear dependence of Whb on metallicity suggested by 
B98 is not confirmed by Figs. |] and ^: there is a break in the 
linear trend at about [Fe/H]~ —1.4. This is more evident for 
thanks to the larger coverage in metallicity of the HST 



B-v 



W KB 
data. 



Figs. | and % shows also that a given change in [Fe/H] 
which is ~ 1.5 times larger than 



B-V 



produces a change in W BB 
that in W^ 1 . 



3. Comparison with theoretical models 

The interpretation of the data was carried out using both the 
VandenBerg et al. (2000) and Girardi et al. (2000) isochrones. 

We first present the two isochrone sets, and then show the 
theoretical prediction for Hhb vs. [Fe/H] together with the em- 
pirical data sets. The thorough inter-comparison is then carried 
out in Sect. f| 



3. 1. The VandenBerg et al. j ^OOQ) isochrones 



The V00 isochrones are computed for 17 values of [Fe/H] 
(-2.3 < [Fe/H] < -0.3), 6 age values (8 < t < 18 Gyr), and 
a fixed mixing length parameter omlt = 1-89. These mod- 
els were computed for a fixed value of [a/Fe] = 0.3, where 
a stands for the alpha-elements O, Ne, Na, Mg, Si, S, Ar, Ca 



Table 1. Measured values of Whb for our sample. An under- 
lined 0.15 dex error in the metallicity, highlights a value of the 
[Fe/H] taken from Harris A 1996b. 
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Fig. 7. V ^ as a function of metallicity, using the HST data 

sample. The CG97 scale is used in the top panel, while ZW84 
is used in the lower one. The dashed lines reproduce the theo- 
retical parameter calculated using the V00 isochrones, between 
8 Gyr (lower line) and 18 Gyr (upper line), and in 2 Gyr steps. 
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Fig. 8. Same as Fig. for the ground-based dataset. 



and Ti. The isochrones have been transformed to the observa- 
tional plane by means of se mi-empirical color-T c fr relations 
(see VandenBerg et al. 2000 ) that are essentially theoretical for 
T e g > 5000 K (where the TD is), and empirically-corrected 
for lower T e g s (where the RGBs are). 



The V00 database includes smooth and well-behaved 
ZAHB sequences for all their metallicity values. They are 
shown in Fig. ||, for both the B — V and the V — I colors. 
All ZAHB sequences have been shifted vertically in this plot, 
so that they coincide at the point where the color is 0, i.e. at the 
turn-down. No color shift has been applied, further confirm- 
ing that the HB-TD point is fixed in color. This figure has also 
been used to define the "mean HB" mentioned previously in 



Sect. 2.1 



The trends with metallicity of W^ v and from 
V00 models (dashed lines) are compared with the data in 
Figs. ^ and [| respectively. As it can be noticed, these mod- 
els reproduce reasonably well the data for W^ v , especially 
if the ZW scale is adopted, but present Wjjjf values that are 
systematically larger than the observed ones by about 0. 1 mag. 
This point is thoroughly discussed in Sect. [4] It is also clear that 
Wb_b increases with age, since the Hayashi track becomes red- 



rV-I 



der as the mass of evolved giants decrease (see also Sect. 4.3) 



3.2. The Girardi et al. (200G) isochrones 



The GOO models were calculated for metallicities Z — 0.0 004, 

0. 001, 0.004, 0.008, 0.019 and 0.03 (Girardi et al. |2000| ). An 
additional set with Z = 0.0001 (or [Fe/H]= -2.28) has 
been computed adopting identical input physics (Girardi, un- 
published). All models have been computed with scaled-solar 
metal ratios and a helium content that mildly increases with Z, 

1. e. Y = 0.23 + 2.75 Z. In this case, the relation [Fe/H] = 
log(Z/0.019) gives [Fe/H] values accurate to within 0.03 dex. 
Thus, the 5 lowest Z values we are going to consider (0.0001, 
0.0004, 0.001, 0.004, and 0.008) would correspond to [Fe/H] 
values of -2.28, -1.68, -1.28, -0.68, and -0.38. The mix- 
ing length parameter was assumed to be q;mlt = 1-68. 

The Girardi et a l. (|200C ) models are computed with OPAL 
(Rogers & Iglesias 1992 ) and Alexander & Ferguson ( 1994 ) 
opacities, and an updated equation of state that considers the 
main non-ideal effects, such as the Debye-Huckel correction. 
In these aspects, their lowest-mass models (oldest isochrones) 
are very similar to many of the non-diffusive stellar model s 
available in the literature (see e.g. Weiss & Schlattl 1998 ). 
However, due to small differences in the adopted physics of 
electron-degenerate matter, and in the way used to generate 
ZAHB structures, Girardi et al. models present ZAHB lumi- 
nosities that are systematically lower (by ~ 0.2 mag) than other 
recent calculations (see Castellani et al. 2000, and Salasnich 



2001, for all details). These difference in luminosity seems to 
be systematic and almost independent of metallicity: in fact, 

GOO ZAHB models obey the relation M^^ yT = (0.163 ± 
0.015) [Fe/H] + constant, (see Salasnich 2001), whose slope is 
consistent with those typically found in most recent models of 
stellar evolution, i.e. 0.18 mag/dex (e.g. Cassisi et al. 1999). 

The GOO core-helium burning models are quite complete 
and well-sampled for masses higher than 0.6 M Q . For lower 
masses, the ZAHB is not as well covered as in V00 models 
and, as a consequence, the identification of the HB-TD is not 
as straightforward. To this aim, we proceeded in the follow- 
ing way: First, we artificially obtained a more detailed ZAHB 
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Fig. 9. ZAHB sequences from GOO models. 



sequence in the HR diagram, by simply interpolating between 
the ZAHB models that have been actually calculated for just a 
few masses M < 0.6 M@. This is displayed in the upper panel 
of Fig. U As it can be noticed, the computed models are dis- 
posed along sequences that are fairly linear in the temperature 
interval logT c ff > 3.8. Thus, we can safely adopt linear inter- 
polations between these different models. Also the mass values 
are linearly interpolated, which represents a less accurate, but 
not critical, approximation. 

Then, the interpolated ZAHB sequences are transformed 
to the observational quantities by adopting the metallicity- 
dependent bolometri c cor rections and color transformations 
from Bertelli et al. ( 1994 ). In the T c g interval we are deal- 
i ng wi th, these transformations are entirely based on Kurucz' 
( 1992 ) library of synthetic spectra. They are mainly a function 
of T c ff and just marginally depend on the surface gravity g. 
Thus, any possible inadequacy in the interpolation of stellar 
masses would not be critical in this step. 

The results for the Afy vs. B — V and My vs. V — I dia- 
grams are shown in the middle and lower panels of Fig. ^ re- 
spectively. Looking at this figure, one can notice that the ZAHB 
sequences of lower metallicities are systematically shifted to 
lower (brighter) magnitudes. 

The position of the turn-down is largely determined by the 
behavior of the V~-band bolometric corrections as a function 
of either T e g or color. We define the turn-down as the point 





[Fe/H] z 



Fig. 10. The W^ v values, as measured from our cluster sam- 
ple (full dots with error bars, and marked with the clusters' 
NGC number), are compared with the values derived from 
Girardi et al. (2000) isochrones (dashed lines). Going from 
above to below, isochrone ages are 16, 14, 12, and 10 Gyr. 

[Fe/H] CG 





Fig. 11. The same as in Fig. but for W^ 1 



of each ZAHB sequence for which B — V = (in the case of 
My vs. B-V diagrams), or V-I = (for My vs. V-I di- 
agrams). Since the zero-points of the "theoretical" photometry 
are such that B-V = V-I = for the AO dwarf Vega, for 
giants of B — V ~ the differences between B — V and V — I 
are of just of a few hundredths of magnitude. Hence, in prac- 
tice, the two different definitions of the turn-down correspond 
to almost (but not exactly) the same ZAHB star. 

Having defined the ZAHB turn-down position for all metal- 
licities, for any isochrone we can measure 

Wife V and ^he 1 
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by simply identifying the color of the point along the RGB that 
is 0.5 mag brighter (in the Wband) than its respective turn- 
down. 



The results for Wjf B v are shown in Fig. 



10 



where 



the theoretical values obtained for isochrones with ages be- 
tween 10 and 16 Gyr are compared to the observational data 
from Sect. |2|. It can be noticed that W BB ~ V values increase 
steadily with metallicity, following the color shift of the RGB. 
If the Carretta & Gratton (1997) metallicity scale is adopted, we 
find a remarkable agreement between models and observations 
throughout the complete metallicity range. The agreement be- 
comes le ss sat isfactory if the data are plotted against the Zinn 
& West (|1984 metallicity scale. 



Similarly, the results for W BB ~ T are shown in Fig. |TT]. 
Again, values are found to increase steadily with metal- 

licity in a way similar to the one present in the data, though (i) 
in this case, there is no data with [Fe/H] > — 1 to compare the 
models with; and (ii) theoretical values are systemati- 

cally larger than the data, by about 0.05 mag . 



4. What determines W H b 

In this section, we will analyze the main parameters affecting 
the theoretical value of Whb- 

The four Figs. ^, [| [l0[ and give us a chance to inter- 
compare two recent sets of theoretical calculations with em- 
pirical data. It is evident from the figures that both sets of 
isochrones reproduce satisfactorily the trend of the observed 
points as a function of [Fe/H] . The agreement is better for the 
BV data. 

For VI, there is a disagreement of ~ 0.05 magnitudes be- 
tween the models and the data, which can be interpreted in two 
ways. Either the models predict a smaller Whb value than ob- 
served, or there is a problem with our measurements. Before 
proceeding further with the discussion, we need to rule out the 
latter possibility. First, we notice that the offset in color is in- 
dependent of the cluster reddening, so it cannot be attributed 
to our differential reddening corrections. Then since Whb is a 
differential quantity, the only justification for such a spurious 
result would be the neglecting of the quadratic (and higher or- 
der) te rms in t he cali bration equations obtained by Rosenberg 
et al. ( |20004 |2000b| ), such that a stretch of 0.05 magnitudes 
over a range of 1 magnitude in color is produced. However, 
the residual scatter around the assumed linear relations was a 
few 0.001 magnitudes over a range in color of more than 1.5 
magnitudes for the standard stars, so we must conclude that the 
problem is indeed in the synthetic colors. 

The main point that emerges from these figures, therefore, 
is that theoretical models seem to describe reasonably well the 
relative position of the HB-TD and RGB, apart from a possible 
zero point difference in V — I. Since the HB-TD color is con- 
stant, Figs. ^, H [l(], and [IT] in fact represent the behavior of the 
RGB color, at about the HB level, as a function of metallicity. 

It is well-known that the RGB color in stellar models de- 
pends on a series of factors. First of all, it depends on metallic- 
ity. Then, it should strongly depend on the efficiency of the en- 
ergy transport in convective envelopes, and hence on the mix- 
ing length parameter a. There must also be a weak dependence 
on the stellar mass, and hence on the adopted isochrone age. 



Some dependence on the model chemistry (helium content and 
degree of a-enhancement) may be present. Finally, the RGB 
color depends on the adopted transformations between T c g and 
color. In the following, we will explore these dependencies in 
order to verify whether the Whb parameter is of some useful- 
ness in setting one of the above parameters. 

4. 1. Dependence of Whb on metallicity 

The dependence on the metallicity of Whb is the most obvious 
and the strongest one. It is clearly visible in Figs. [}] and |j. As 
the metallicity increases, the RGB becomes redder and redder, 
while we do not expect a variation in the HB-TD color, and 
therefore Whb increases. The effect is stronger in (B — V) than 
in (V — I). We will discuss this dependence in more details in 
Sect. 6. 

4.2. Dependence of Whb on the mixing length 
parameter a 

Among the various parameters fixing the RGB position, the 
mixing length parameter a plays one of the main roles, since 
in RGB stars a significant fraction of the energy flux is trans- 
ported by convection. The larger the value of a, the higher the 
temperature of the RGB. 

In order to evaluate the effect of a in the Whb parameter, 
we computed an additional set of stellar models, from the zero 
age main sequence (ZAMS) up to the He-flash, for a equal to 
1.30 and 2.00, respectively. Only stars of 0.8 A/© were consid- 
ered; this mass value is close the one found in the upper RGB 
of isochrones that are 10 to 15 Gyr old. Together with the set of 
a = 1.68 tracks already available (GOO), the new tracks allow 



us to evaluate how W BB v and W BB 1 change as a function 
of a. The W BB V and W BB ^ values for all our models with 
0.8 Mq, are plotted in Figs. [12] and [l3|, together with the values 
derived from our 10 — 16 Gyr old isochrones. 

It can be noticed that, for lower metallicities, the points 
corresponding to 0.8 Mq, a = 1.68 models, coincide with 
the locus of the isochrones. For higher metallicities, this does 
not happen, since metal-rich 10 — 16 Gyr isochrones typically 
present RGB stars of higher masses (0.9 — 1.0 Mq), and hence 
hotter and bluer than 0.8 Mq RGB stars. 

Other points to be noticed are that: (i) Models calculated 
with different a present virtually the same value of core mass 
at the helium flash as the a = 1.68 ones, and a similar amount 
of dredged-up helium. This implies that the ZAHB models de- 
rived from these tracks would present essentially the same lu- 
minosities, (ii) Stars hotter than T e g — 7000 K have their radii 
essentially insensitive to the treatment of the outer convection, 



and therefore independent of a (see Renzini & Fusi Pecci 1988 



and figure 1 in Castellani et al. g_999(). Hence, even if we did not 
compute ZAHB models with different a values, we know that 
the HB-TD position would not have changed. 

Then, we can use the computed 0.8 Mq RGB models 
to derive the dependence of Whb on a: For each metallic- 
ity, we have fitted straight lines to the W BB V versus a, and 
Wjib" 7 versus a data, using the 3 RGB tracks we have for 
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Fig. 12. Open squares: W^ v values measured for a set of 
0.8 Mq evolutionary tracks of varying metallicity and mixing 
length parameter a. For each metallicity, a values are 1.30, 
1.68, and 2.00, going from above to below. For the sake of 
comparison, we also plot the W^ v values as derived from 
the same isochrones as in Fig. [l(^ (dashed lines), and the obser- 
vational values for our cluster sample (small dots). 



each metallicity. In all cases, the linear fit produced an excel- 
lent description of the data. Therefore, we can conclude that the 
1 (a) relations are very much linear. The 



w£ B v (a),<mdW" B 



slopes of the fitted lines, AWjf~ v / Aa and AW^ 1 /Aa, are 
presented in Table ||, for each [Fe/H]. 

The numbers in Table ^ confirm the strong dependence of 
Whs on a, to be compared with the much lower dependence 
on age. 



4.3. Dependence ofWnB on age 

The dependence of Whb on age is already illustrated in Figs. ^, 
H [To| and [jr], where isocrohones for a wide range of ages are 
presented. In Table ^, instead, we tabulate the derivatives of 
Whs with respect to age, at 14 Gyr, as derived from GOO mod- 
els. Surprisingly, in V00 models these derivatives are about 
twice as large; but anyway, it is clear that M^hb changes very 
little with age. Taking int o acco unt the small age dispersion 
found by Rosenberg et al. ( 1999 ) for the same clusters used in 
the present investigation, we conclude that the effects of age 
variations from cluster to cluster are much smaller than the er- 
ror bars of the single data points. 



4.4. Dependence ofWnB on the helium content 

When the helium content is increased, the HB luminosity also 
increases according to the law AV/AY = —3.22 for 0.23 < 




Fig. 13. The same as in Fig. H], but for W^ 1 . 



[Fe/H] Cl 




Fig. 14. The dashed lines show the effect, on W^ v ', of vary- 
ing the helium content on the theoretical models. See text for 
the explanation. 



Y < 0.27, i.e. the values suggested by current observations (R 
parameter, HII regions, etc.). This means that, when measuring 
Whb, we will consider a brighter RGB point as well. This point 
will be redder, whereas the HB-TD does not change in color, so 
the net result is that increasing Y will lead to a larger W-rb- The 
size of this effect was measured using a 12 Gyr isochrone from 
V00, varying [Fe/H] in the whole range, and for five values of 
Y. 

The comparison with the observations is shown in Figs. [IJ 
and O, for the HST and ground-based samples, respectively. 
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Table 2. Derivatives of Whb cf. GOO models. 
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Fig. 15. Same as Fig. |l4|, for the ground-based observations. 

The figures show that, while there is a zero-point problem 
(which could be solved changing the reference age), the rela- 
tive trend is well reproduced by the models. In any case, vari- 
ations in helium content have a negligible effect on VFhb ■ The 
dispersion of the empirical points due to the observational un- 
certainties is larger than the dispersion we would expect from 
a reasonable cluster to cluster variations in the helium content. 
Therefore, we conclude that the dependence of M'hb on the He 
mass fraction can be ignored. 

4.5. Dependence of W^b on a-enhanced metal 
ratios 

Halo populations are expected to present a-enhanced metal ra- 
tios, with [a/Fe] ~ 0.3 dex. For a given metal fraction Z, en- 
enhanced models have an iron content [Fe/H] that is depleted 
by a factor that corresponds roughly to the degree of enhance- 
ment [a/Fe]. Moreover, for sufficiently low metallicities, a- 
enhanced evolutionary tracks can be replaced by th eir scaled- 
solar counte rparts of same Z (see Salaris et al. 1996| ; Salaris & 
Weiss |1998| ). 

Therefore, for low metallicities we can expect that mod- 
els with solar-scaled abundances could be used to reproduce 
a-enhanced isochrones, provided that their [Fe/H] values are 
changed by about —0.3 dex. Of course, this is only a first- 



order approximation to the question, especially for the models 
of hig her metallicities ([Fe/H] ^ 0.004, see Salasnich et al. 
2000 ). Another approximation is in the fact that the T c ff-color 
transformations in use have been derived from model atmo- 
spheres of scaled-solar composition, and not from a-enhanced 
ones (which are not yet available). Whether this is critical, is 
yet to be investigated. 

Keeping these points in mind, we can now get some in- 
sight on how much Whb depends on the a-enhancement, by 
taking advantage of the fact that V00 isochrones are enhanced 
while those of GOO are not. As a caveat, we must stress that 
ideally one would like to measure the size of the effect on oth- 
erwise identical models, since e.g. different color transforma- 
tions could mask the true trend. As a first order approach, we 
will nevertheless try to understand what happens if we inter- 
pret the differences in the two isochrone sets purely in terms of 
a-enhancement (i.e. as just an offset along the abscissae). 

Let us then compare Fig. @ to Fig. [fi| and H to looking 
for example at their [Fe/H] position at Wub = 1- We can see 
that a-enhanced isochrones look indeed more iron poor than 
solar-scaled isochrones. In particular, the displacement for the 
(B - V) color is A [Fe/H] ~ -0.1, while it is A [Fe/H] ~ 
—0.2 for the (V — I) color. The effect then goes in the right 
direction, although the size of the offset in the iron scale is less 
than expected (but remember the above caveat). 

If we now introduce the data, we see that the (V — I) color 
is the one that gives more troubles. There is a certain degree of 
agreement if we consider the GOO theoretical trend, when data 
are plotted on the ZW84 metallicity scale. In all the other cases 
the isochrones look too iron poor, regardless of the metallicity 
scale. 

The case of the (B — V) color is less clear-cut, since the 
choice of the metallicity scale now plays an evident role. If we 
assume that the CG97 scale is correct, then the GOO models 
show the best agreement, so it looks like no a-enhancement 
is required. On the contrary, the V00 models are the ones that 
better reproduce the observed trend on the ZW84 metallicity 
scale. The GOO models should then be enhanced in order to 
reach the same agreement. 

In conclusion, both theoretical calculations show problems 
in reproducing the (V — I) trend of M'hb vs. [Fe/H], most 
probably due to the color- T c ff relations (see also Sect. 4.6). 
If we then just rely on the (B — V) color, then it is really a 
matter of choosing one's preferred metallicity scale. Since a 
elements are indeed enhanced in GGCs, and trusting the mod- 
els, one should choose the ZW84 scale and V00 isochrones. 
It is also likely that GOO isochrones will be consistent with the 
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data on the ZW84 scale, once the a-enhancement is introduced. 
However, it looks like the real way out should be a third inde- 
pendent determination of the GGC metallicity scale. 

If the ZW84 scale could be demonstrated as the correct one, 
then there would not be any need for a revision of the present 
picture (besides some revision of the (V — I) -f- T c g relations). 
If instead the CG97 scale is the good one, and retaining the a- 
enhancement scenario, then even (B — V) color-temperature 
relations should be revised. 



4.6. Dependence ofWnB on the T cS -color 
transformations 

We have already noticed that the models that so well fit the 
Wife v data (Fig. |l0|), do not fit as well the W^ B 1 (Fig. [ll|) 
ones. In particular, while the theoretical trend seems to be the 
same as the observed one, there is a zero point shift for W^ 1 . 
Since the isochrones we are using in these plots are the same, 
this is probably indicating that we have problems with the color 
transformations, either for B — V or V — I, or both. In fact, dif- 
ferences between Kurucz and empirical T c s vs. color transfor- 
mations have been noti ced by several authors ( e.g. G ratton et 



al. 1996; Castelli et al. 1997; Weiss & Salaris 1999). Notice 



that, in our isochrones, it would be enough to have the RGB 
just 0.05 mag bluer in V—I to find a perfect agreement be- 
tween models and observations. A color shift of this order can 
be caused even by the different filter transmission curves used 
by several authors. 

Just to give an order of magnitude to the uncertainty due to 
the color transformations, in Table || we present the changes in 
color, A(B — V^rgb and A(V — /)rgb, that we obtain for 
the RGB point 0.5 mag above the TD level in GOO 12 Gyr-old 
isochrones, in the case we adopt Al onso e t al. ( 1999 ) T c $ -color 
transformations instead of Kurucz ( 1992 ) ones. Unfortunately, 
such a difference can not be found for all values of [Fe/H], sim- 
ply because at low met allicit ies we go out of the applicability 
range of Alonso et al. ( 1999 ) formulas (see their tables 2 and 
3). 

As can be noticed, for higher metallicities we find consis- 
tent shifts (~ 0.15 mag) in the B — V color of the RGB; how- 
ever, it is not possible to establish how these differences depend 
on metallicity, since only 2 values were derived (for Z = 0.004 
and 0.008). For V — I, instead, the differences are surprisingly 
small (~ —0.04 mag), and seem not to depend on metallicity. 

This exercise indicates that the uncertainties in T e ff-color 
transformations could sensibly affect our results for W^ v , 
but not for . However, we remark that a comparison is 

essentially missing for low metallicities. 

We note that Alonso et al. ( 1999) ) constitutes the most up- 
dated source of empirical T e ff-color relations for red giants. 
The comparison with other different transformations from the 
literature would not help much to clarify the issue of the behav- 
ior of Whb- 



5. Can Wub be used to calibrate the mixing 
length parameter a ? 

In the previous section we showed that the mixing length pa- 
rameter a plays a major role in determining the values of Whb, 
being by far more influent than the cluster age and helium con- 
tent. 

Moreover, a rapid inspection of Figs. |l2| and [jj] suggests 
that present data are confined within a relatively narrow range 
of a values. A natural question that comes to us, then, is: Is 
the data compatible with a single value of the mixing length 
parameter a ? And, if yes, can we then use the data to single 
out a "best value" for a ? 

In answering these questions, however, we should keep in 
mind the uncertainties in the T e ff-color transformations, and 
the enhancement of a-elements, that are the main uncertainties 
in the comparison between models and data. They have been 
shortly mentioned in the previous section. Regarding the color 
transformations, we notice a mismatch in V—I data (Fig. pj] ), 
and the possibility that similar effects are present also in the 
B — V color. Considering this, we conclude that a single "best 
fitting value" for a cannot be derived, unless we know the color 
transformations with an accuracy of some hundredths of mag- 
nitude. The best we can do, for the moment, is to check whether 
present models deviate from the assumption of a single a being 
valid throughout the entire [Fe/H] range of our observations. 
In this respect, the most important conclusion of Sect. \.d is 
that the offset introduced by changing the color transformations 
seems not to depend on metallicity. 

In order to make this check, we first measured the distance 
of each data data point from a given set of isochrones (for 
a fixed age and a), AWhb- Then, we translated the AWhb 
values into differences in alpha, Aa, using the AWnB/Aa 
derivative appropriate to each metallicity (Table |2|). Then, any 
trend of Aa with metallicity can be interpreted in terms of 
changes in the mixing length parameter, under the assumption 
that no other parameter is playing a significant role, as we have 
demonstrated in the previous section. 

This was performed separately for AW^ V and AW^ 1 . 
Figs. [If] and [[7] show the Aa values as a function of metal- 
licity, obtained having, as the reference, the GOO isochrones of 
age 14 Gyr. Noticeably enough, the differences in a show no 
significant trend with metallicity. The only case in which some 
marginal trend of this kind could be present, is the one for V—I 
data on the ZW scale, for which AW^ 1 seems to slightly de- 
crease for [Fe/H] > —1.5. 

We have then modeled the data in Figs. |l6| and |l7| by least- 
squares fitting of either (i) a constant value, or (ii) a straight 
line. The results are presented in Table ||. 

In the case of BV data, the fitted lines turned out to have 
a negligible slope, and there is virtually no improvement in the 
fittings as we pass from a constant to a straight line. The same 
applies for the VI data in the CG97 scale. The case of VI data 
on the ZW scale, with its mild slope (— 0.17±0.07), constitutes 
the only exception. 

Perhaps more importantly is the fact that, assuming that a 
is constant (first part of Table ||), we get extremely small values 
for its dispersion, i.e. a a ~ 0.023. This has to be compared to 
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Table 3. Fitting parameters for the data in Figs. ^6] and |l7| 



(i) A constant value: 


Aa = a 










Case 


a 








x 2 


CG scale, BV data 


-0.036 


0.022 






10.28 


7W scale BV data 


-0.157 


0.025 






9.82 


CG scale, VI data 


0.178 


0.022 






13.11 


ZW scale, VI data 


0.133 


0.024 






15.59 


(i) A straight line: A 


a = a + b ([Fe/H] + 1.5) 






Case 


a 




b 


Ob 


x 2 


CG scale, BV data 


-0.036 


0.023 


-0.009 


0.059 


10.26 


ZW scale, BV data 


-0.157 


0.025 


-0.027 


0.064 


9.64 


CG scale, VI data 


0.174 


0.024 


-0.028 


0.066 


12.92 


ZW scale, VI data 


0.108 


0.026 


-0.173 


0.071 


9.71 



Lrc/iij c 




[Fe/H] CG 



Fig. 16. The estimated Aa differences for each object observed 
in B — V, and the reference GOO isochrone of age 14 Gyr, as a 
function of metallicity in the ZW (bottom panel) and CG97 
scales (upper panel). Notice the absence of any significant 
trend with metallicity. 



typical values of a assumed in evolutionary calculations, i.e. 
a ~ 1.7. Overall, these numbers indicate that, indeed, the em- 
pirical data can be very well approximated by a constant value 
of a. 



6. A new metallicity index? 

We have shown that there is a relation between VFhb and the 
metal content of the cluster. Since Whb has only a second or- 
der dependence on the reddening (which can be corrected as 
described in Appendix), it is a potentially interesting metallic- 
ity index. We have verified whether it is possible to use this 
parameter as a new metallicity index. In Fig. [l^ we fitted a 



< 0.3 





Fig. 17. The same as in Fig. [T(| but for the clusters observed in 
V — I. 



parabola to the observed vs [Fe/H] relation for the 

HST sample, obtaining a 0.1 dex rms residual. 

The sensitivity of the parameter, although lower than that of 
other traditional indices, can be used to obtain a first estimate 
of the metallicity. A typical 0.04 magnitude error on Wjjb" 
would translate into a ~ 0.2 dex uncertainty on [Fe/H]. On 
the other hand, the lack of metal-rich clusters, and the much 
shallower dependence on metallicity of this parameter, prevents 



a reliable calibration of W, 



v-i 

HB 



There could be some concern about the inclusion of the two 
metal-rich clusters (NGC 6388 and NGC 6441) in the B - V 
calibration, since it is still not known what drives the peculiar 
morphology of the HB of these two objects. Therefore, we have 
also calibrated the W^ v index in the -2.1 <[Fe/H]< -1.0 
metallicity interval. In this case we used the weighted linear 
least squares fit, which is shown in Figs. |l^ and The result 
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[Fe/H] =a-W ? +b-W+c 



-0 5- [ pe / H ] = o ■ w 2 + b • w 



o = -6.20 
b = 14-98 
c = -9.52 
rms = 0.10 





Fig. 18. Calibration of W® B V as a metallicity index. The ob- 
served data points have been interpolated with a second order 
polynomial. 





Fig. 19. Calibration of W^ v for the HST clusters, obtained 
with the exclusion of NGC 6388 and NGC 6441. The param- 
eters of the weighted least square interpolation are reported in 
the figures, both for the CG97 (upper panel) and ZW84 (lower 
panel) metallicity scales. 



of the fits shows that the relation for W^ v nas still a for- 
mally low dispersion (0.1 dex rms), whereas a large 0.3 dex 
dispersion is shown by the residuals for the parameter. 
Again, this confirms the higher sensitivity of W HE j~ to metal- 
licity changes. 

In Fig. |A.l[ which shows the trend of the second order red- 
dening corrections that must be applied to Whb> the right verti- 
cal axis shows how a difference in Whb translates into a differ- 



Fig. 20. Same as Fig. hjJ, for the ground-based data. 



ence in metallicity, using the above relations. It clearly shows 
that, especially for the highly reddened clusters, neglecting the 
differential reddening corrections would lead to significant sys- 
tematic errors in the estimated metallicity. For instance, the er- 
ror on [Fe/H] for NGC 6287 would be of the order of 0.2 dex. 
Naturally, these errors are larger for metallicities obtained us- 
ing the V — I calibration. 

In any case, we must note that: (i) The rms error are just 
formal fitting errors, while the real uncertainty on the metallic- 
ity is higher, as it must take into account the error on the metal 
content of the calibrating clusters and the photometric calibra- 
tion errors; (ii) The relations obtained in this paragraph can be 
used only in a limited metal interval and for a limited sample 
of clusters, which must have a blue HB. 



7. Summary and discussion 

In this paper, we took advantage of the homogeneous pho- 
tometric databa ses of Galactic globular cluster s, pre sented in 
Rosenberg et al ( |2000a[ |2000b| ) and Piotto et al. ( |2002| ), to com- 
plete another step in the characterization of the morphological 
properties of the CMD of GGCs, and in the fine-tuning of the 
theoretical models. To this aim, the Whb parameter, originally 
defined in B98, has been measured for all suitable clusters, and 
appropriate corrections to account for differential reddening ef- 
fects have been applied. The quality of the new data sets has al- 
lowed to settle the original question posed in B98. The trend of 
Whb with metallicity showed a dispersion that was larger than 
formal measurement errors, a fact that can now be attributed to 
the inhomogeneous data sources employed. 

The dispersion of Whb around the mean trend with metal- 
licity, is now compatible with the error bar on the data points. 
This means that (a) one can reverse the argument and use Whb 
to have a first-order guess on a cluster metallicity (Sect. ||), and 
(b) that whatever other variables influence the parameter, they 
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must have a well-defined dependence on [Fe/H] (including zero 
dependence). 

The second point was investigated in some detail by com- 
paring the data to two independent sets of theoretical calcula- 
tions. First, it was noted that the observed trend of Whb with 
metallicity (which is stronger for the B — V color) is well re- 
produced by both isochrone sets, although a zero-point prob- 
lem in the (V — I) color-temperature transformations seems to 
be present. 

In order to test the dependence on all other variables influ- 
encing Whb, their values were varied within reasonable limits. 
For any variable x, we checked the Whb vs. x trend at several 
fixed metallicities, and noticed that these trends can be well 
approximated with linear relations. This means that we could 
use the slope AWhb/Ax to rank the relative dependencies. 
We concluded that the influence of the helium content is neg- 
ligible, and that that of the mixing-length parameter a is much 
stronger than that of the age (see Table ||). This fact was then 
exploited to investigate a long-standing theoretical problem of 
stellar evolution, whether a should be varied with metallicity 
or not. From our comparisons, we find no trend of a with the 
metal content of a cluster. 

With respect to the color- T c s transformations, a test was 
made by adopting either the Alonso et al. (1999) or the Kurucz 
( 1992 ) relations. This showed that, potentially, the B~V trans- 
formations could make a much larger difference than the V — I 
transformations. However, since actually the theoretical V — I 
colors show the greater problems, the agreement between the 
two transformations implies that a deep revision for them is 
needed. 

Finally, we examined the question of whether models with 
enhanced a-elements better reproduce the CMD morphology 
of GGCs. The conclusion is that, unfortunately, the existence 
of two discrepant metallicity scales leaves this question open. 
The best results are obtained in the a-enhancement scenario, 
and adopting the ZW84 metallicity scale. However, one could 
well adopt the CG97 scale and claim that the B — V color- 
transformations are wrong, since we have seen how much they 
change from author to author. Thus, we must conclude that 
there is an urgent need for an independent metallicity scale 
(which also involves the question of g alacti c chemical evolu- 
tion models, see Saviane & Rosenberg 1999). 
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Appendix A: Differential reddening corrections 

A. 1. Corrections according to Olson 



Olson (1975) obtained the following relation: 
R = 3.25 + 0.25 x (B - V)o + 0.05 x E B -v 



>From the definition of R = 
we can write: 



and of the color excess, 



W? 



HB 



(B - V) RGB — (B — V) 



o 

TD 



= (B- V) RGB - E 



RGB 
B-V 



obs 



-{B 


- V)td + 


&B-V 


"obs + 


( fTD pRGB \ 
\&B-V ~ &B-V) 


(A v 


_ Av 




V R 


TD R 


rgb) 



Since Av is little dependent on the color, we substitute an av- 
erage value in the equation, (Av), obtaining: 



HB 



obs 



(Av) 



R 



RGB 



R 



VD 



R TD . R RGB 



Now we substitute the product of i? TD times i? RGB by R 2 , and 
using Olson expression, we find: 



W bs+ < A v > 



0.25 [(B - lOo - (B - VOo ] 



<b 

TDl 



= W ohs 



R 2 



0.25E B -yA{B - V) 
R 



Here Eb-v f = ^ ~ ) is the average color excess of the single 
cluster, taken from Harris (119961 ), and we assumed R = 3.1. As 
a typical value of A(B — V)o we used 0.8 magnitudes, which 



is the W£n V value of NGC 1904, taken from B98 who ap- 
plied the same differential reddening correction. In conclusion, 
the observed values of must be corrected with the ex- 

pression: 

W£ B = W ohs 



0.064 x Eb-x 



(A.l) 



A.2. Corrections according to Grebel and Roberts 

Grebel and Roberts ( |1995| ; GR95) provide tables that give 
color excesses, absorptions, and the R ratio for different pho- 
tospheric temperatures, gravities and metallicities ([Fe/H] = 
0.0, —1.0, —2.0). Since we are applying a second-order cor- 
rection, we used the tables relative to the intermediate metal- 
licity. 

We start again from the relation 



HB 



obs 



(E 



TD 
B-V 



pRGB 
^B-V 



And we find the color excesses using the relation 

< A v > 



R = 



/(tab 
A V 



E 



'tab 
B-V 



Z?col 
^B-V 



(A.l) 



(A3) 



where Ay h and are the tabular values from GR95, 

E^iy is the color excess for a given cluster and given color 
(to be determined), and (Av) is an average value of the total 
absorption, which can be calculated with the expression: 



< A 



>= 3.1 x E B at _ v 



where E^Iy is the mean reddening taken from Harris ( |1996[ ) 



Equation (A. 3 ) is justified by the fact that the ratio of total 
to selective absorption is constant. Thus, it can be used to find 
Eb-v for any TD and RGB color, since: 



peat 

^B-V - 3-1 R — 



(A.4) 
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We now insert Eq. (A.4) into Eq. (A.2), and taking R from 
GR95 tables, we obtain 



W^HB 



W obs + 3.1 x 



C X E^y 



1 



1 



R 



I'D 



R 



RGB 



peat 
x &B-V 



(A.5) 



This last equation (A.5) is formally similar to Eq. (A.l ), but in 
this case c is a variable that depends on the color of the two 
CMD points. The variations of c are, as expected, of the order 
of a thousandth of a magnitude. 

A3. Corrections to the (V- 1) measurements 

The analogous of Eq. (|A.2|) is the following 



W£ B = W obs + (E 



TD 

V-I 



E 



RGB-, 
V-I ) 



The Ev-i values can be computed from Eb-v 



E 



col 
V-I 



(E v - 



\E B - 



t ab 



x E 



B-V 



(A.6) 



(A.7) 



where, like before, "tab" quantities depend on the temperature 
and are taken from GR95 tables, while "col" quantities are re- 
ferred to different zones of each cluster's CMD. Thus E 



the value of the blue color excess (Eq. |A.4| ), while Ey _j is the 
value of the color exc ess t hat must b e co mputed. 

Substituting Eq. (A.7) into Eq. ( |A.6| ), and using Eq. (A.4). 
one obtains: 



col 
B-V 



w ob8 



( Ey-i 
\E B -v 



tab 



^B-V 



TD 



w, 



obs 



3.1 



1 E V - 



R Eb- 



B-V/td 



tab 



Ey-i 

E B -v 
1 E v ^j 



tab 



<B 



ttiRGB 
&B-V 



R Eb- 



RGB 
tab 



B-V J RGB 



77'Cat 

&B-V 



w 



bs 



77'Cat 
^B-V 



The corrections to the two Whb are plotted in Fig. A.l 



as a 
- 2 



function of Eb-v- It is then clear that (V — I) colors are 
times less sensitive to this effect. 
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